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Abstract 

We propose a new blow-up criterion for the 3D Euler equations of incom¬ 
pressible fluid flows, based on the 3D Euler-Voigt inviscid regularization. 
This criterion is similar in character to a criterion proposed in a previous 
work by the authors, but it is stronger, and better adapted for computa¬ 
tional tests. The 3D Euler-Voigt equations enjoy global wcll-posedness, and 
moreover are more tractable to simulate than the 3D Euler equations. A 
major advantage of these new criteria is that one only needs to simulate the 
3D Euler-Voigt, and not the 3D Euler equations, to test the blow-up criteria, 
for the 3D Euler equations, computationally. 
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1. Introduction 

A major difficulty in the computational search for blow-up of the 3D in¬ 
compressible Euler equations is that one must seemingly simulate the 3D 
Euler equations themselves to obtain information about singularities. Near 
the time of a potential singularity, sufficient accuracy of simulations of the 3D 
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Euler equations can be challenging to obtain, due to the need to resolve spa¬ 
tial derivatives that are potentially infinite. However, two new blow-up cri¬ 
teria, one proved by the authors in [[32], and another provided in the present 
work, provide a path around this difficulty by using the Euler-Voigt inviscid 
regularization. In particular, by tracking the L 2 -norm of the vorticity of so¬ 
lutions to the 3D Euler-Voigt equations, as a certain regularizing parameter 
a tends to zero, these new criteria allow one to gather evidence for potential 
singularities of the 3D Euler equations by only simulating the 3D Euler-Voigt 
equations. This is advantageous from a computational standpoint, because 
the L 2 -norm of the spatial gradient of solutions to the Euler-Voigt equations 
is uniformly bounded in time, for any fixed value of the regularization pa¬ 
rameter a. Furthermore, all higher-order norms grow at most algebraically 
in time [32|, which implies that pointwise spatial derivatives grow at most 
algebraically in time. No such results are known for the 3D Euler equations. 
This means that simulating the 3D Euler-Voigt equations is computation¬ 
ally more tractable than simulating the 3D Euler equations, since achieving 
sufficient accuracy requires that simulations have high enough resolution to 
resolve spatial gradients. 

In this work, we provide a new blow-up criterion that is similar in charac¬ 
ter to the criterion in (32j, but that has several advantages over the previous 
criterion. Our main focus is on differences in the computational implemen¬ 
tation of the two criteria. However, one analytical advantage is that the new 
criterion is potentially stronger than the previous criterion. This is because 
the set of singularities it can detect is a (possibly proper) superset of the 
singularities detectable by the criterion in [32]. From the standpoint of com¬ 
putational implementation, the criterion in |32| does not allow for simulations 
with adaptive or variable time-stepping, requiring a fixed time step in each 
simulation, and also agreement in time-steps across all simulations as the 
regularization parameter varies. (Interpolating in time is also an option, but 
introduces additional approximation.) The new criterion only requires that 
the simulations end at a common final time. Furthermore, the previous crite¬ 
rion requires computation and output of the L 2 -norm of the velocity gradient 
(or the vorticity) at every time step, which can be costly, requiring additional 
memory storage, and—in parallel simulations—additional communications. 
The new criterion requires this data only at the final time. We note that 


1 Of course, it may be that there are no singularities in the 3D Euler equations. 
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both of these criteria are only known to be sufficient for blow-up; they are 
not known to be necessary for blow-up, unlike, e.g., the Beale-Kato-Majda 
criterion [2|. Currently, both criteria are being tested computationally plj ]. 
The computational search for blow-up of the 3D incompressible Euler 

, 23, 0], and the refer- 


equations has a long history (see, e.g., [14, |19|, 21 


ences therein). Traditionally, one attempts to identify singularities by means 
of blow-up criteria based on quantities arising from the 3D Euler eq uations. 
There are many such criteria in the literature (see, e.g., [2, 10, EH, 00, S3, 
and the references therein). Computational tests of the type of blow-up crite¬ 
ria described here and in [32] require more simulations (of the 3D Euler-Voigt 
equations) than tests requiring a single simulation of the 3D Euler equations, 
since one must run simulations for several values of the regularizing parame¬ 
ter a. However, these simulations require less resolution than simulations of 
the Euler equations as discussed above. 

The Euler-Voigt inviscid regularization of the Euler equations is given by 


— a 2 d t V 2 u + d t u + (u ■ V)u + Vp = 0, 

V • u = 0, 
u(x, 0) = Uo(x). 


(l.la) 

(lib) 

(1-lc) 


The parameter a > 0, having units of length, is the regularizing parame¬ 
ter. Formally, setting a = 0, we recover the incompressible Euler equations. 
The fluid velocity field, u = u(x, t), and the fluid pressure, p = p(x, t) are 
the unknown quantities. We consider system (II.ip in a periodic box D := 
M 3 /Z 3 == [0, l] 3 . We assume that the spatial average of f Q u 0 (x) dx = 0. 
With (11.lap , this implies f n u (x,t) dx = 0 for all t. From now on, we denote 
by u“ the solution to (II.ip . and by u a solution to the Euler equations (i.e., 
(11.11) with a = 0), both starting from the same sufficiently smooth initial 
condition u 0 . We denote the vorticity w := V x u, and uj a := V x u". 

The Euler-Voigt equations were first described and analyzied in [E], where 
they were shown to be globally well-posed for all initial data u 0 G H 1 and 
all t > 0. We note that their viscous counterpart, called the Navier-Stokes- 
Voigt equations, were proposed and studied much earlier in [38, [33], as a 
model for Kelvin-Voigt viscoelastic fluids. The Euler-Voigt equations have 


been studied computationally in [15|,|31]- The Euler-Voigt and Navier-Stokes- 


Voigt equations, along with extensions of these models, have been studied in 
both a nalytic al and numerical contexts (see, e.g., 0 0 B 0 0, E3 0, 0, 


0,0, H 0 0,0,0,0, EH). 
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The following theorem was proved in 32j (see also a similar theorem for 
the surface quasi-geostrophic (SQG) equations in [28]). 


Theorem 1.1 ([320. Assume u 0 e H s , for some s > 3, with V ■ u 0 = 0. 
Suppose there exists aT* > 0 such that the solutions u“ of (II.If) . with initial 
data u 0; satisfy 


sup limsup (o;||Vu“(f)|| L 2 ) > 0. (1.2) 

t€[0,T*] a->0+ 

Then the 3D Euler equations, with initial data u 0; develop a singularity within 
the interval [0, T*]. 

Remark 1.2. Since V ■ u“ = 0, integration by parts can be used to show that 

||Vu Q (t)|U 2 = ||a;“(t)|| L2 . (1.3) 

Therefore, in line with other blow-up criteria in the literature, (j 1.2 j) can be 
seen as a condition on the vorticity, albeit from the 3D Euler-Voigt equations 
rather than the 3D Euler equations. 

A technical difficulty arises in computational tests of Theorem ll.il Math¬ 
ematically, one may imagine fixing a t > 0 and computing 

limsup (a||Vu a (f)||i, 2 ). (1.4) 

a— » 0 + 

However, computationally, it is more natural to first fix an a > 0, as a pa¬ 
rameter, and then to compute u a (t) as t increases up to a time T (e.g., by 
a standard time-stepping method). Therefore, to construct curves of a vs. 
a|| Vu“(f)|| L2 for each fixed t, one must jump from solution to solution as a 
varies. This gives rise to some of the technical issues discussed above. How¬ 
ever, suppose for a moment that one is allowed to commute the two limiting 
operations in (11.21) . In this case, one would then be interested whether 

limsup (a sup ||Vu“(f)|| L2 N ) > 0. (1.5) 

a—>0+ ' te[0,T*] 2 

The quantity in (jl.5[) is arguably easier to track, as discussed above. It is 
the purpose of this work to show rigorously that (II.5ft implies that the 3D 
Euler equations develop a singularity within the interval [0,T*]. 
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2. Notation and Preliminary Results 


We denote by L p and H s the usual Lebesgue and Sobolev spaces over the 
periodic domain Q = [0, l] 3 := M 3 /Z 3 , respectively. It is a classical result 
(see, e.g., 01, 35[) that, for initial data u 0 E H 3 satisfying V • u 0 = 0, a 
unique strong solution u of the 3D Euler equations exists and is unique on a 
maximal time interval that we denote by [0,T*). Moreover, one has 


|u(t )|| L 2 = ||u 0 ||l 2 on [0, T*). 


( 2 . 1 ) 


Equation (12.11) holds under weaker conditions on the smoothness of the solu¬ 
tions of the 3D Euler equations, as it was conjectured by Onsager (see, e.g., 
0 , However, the existence of such weak solutions for arbitrary 

admissible initial data is still out of reach. In [lj , it was shown that a certain 
class of shear flows are weak solutions in L°°((0, T); L 2 ) that conserve energy. 
Furthermore, families of weak solutions that do not satisfy the regularity as¬ 
sumed in the Onsager conjecture have been constructed that do not satisfy 
(12.21) QMM 2^1 -The following “o-energy equality” was proven in Q. 


Theorem 2.1. Let uo E H 1 with V-uo = 0, and let u“ be the corresponding 
solution to CCD. Then, for any t E M, 


u 


£2 + a 2 ||Vu c 


2 11 112 1 211V7 112 

L 2 — 11 U-0 11 z.2 + OL || VUo|| L 2. 


( 2 . 2 ) 


The following convergence theorem was proven in 32 


Theorem 2.2. Let u 0 E H s , s > 3 with V ■ u 0 = 0, and let [0,T*) be the 
corresponding maximal interval of existence and uniqueness of the solution, 
u, to the 3D Euler equations. Choose T E [0,T*). Then there exists a 
constant C > 0, which depends on sup 0<t<T ||-u(t)||^3 7 such that for all t E 

[o,U 


||u(() - + a 2 ||V(u(«) - u“(t))|||j < C o 2 (e°' - 1). (2.3) 

3. An Improved Blow-up Criterion 

Let T > 0 be given. Assume that a given solution to the Euler equations 
is smooth on [0, T], so that in particular, (12.11) holds. We emphasize that 
(12. lj) depends on the regularity of the Euler equations, and if a finite-time 
singularity develops, (12. ip might not hold. 
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Theorem 3.1. Let u 0 G H s , s > 3, with V • u 0 = 0, and let u“ be the 
corresponding unique solution of (11.11) . Suppose that 

limsup sup a||V u“(t)|U 2 >0, (3.1) 

a—>-0+ te[0,T] 

for some T > 0. Then the unique solution to the 3D Euler equations, with 
initial data Uo, must develop a singularity within the interval [0, T]. 

Proof. We prove the contrapositive. Assume that u is a solution of the 3D 
Euler equations, with initial data Uo G H s , s > 3, that remains smooth on 
the interval [0,T]. In particular, the smoothness implies that (12. ip holds. 
From (12.31) there exists a constant C > 0, depending on sup 0<t<T ||«(t)||#3, 
such that 

||u“(t )|| z ,2 > ||u(f )|| L 2 - Ca(e ct - 1) 1/2 > ||u(f )|| L 2 - Ca(e CT - 1) 1/2 (3.2) 
= ||u 0 || L 2: - Ca(e CT - 1) 1/2 . 

Here, we have used (12.11) . Let a > 0 be small enough that the right-hand 
side is positive (i.e., a < ||u 0 || L 2 /(C(e CT — 1) 1//2 ). Squaring, we obtain 

||u“(t )||^2 > Huollia - 2C'a||u 0 || L2 (e c,T - 1) 1/2 + C 2 a 2 (e CT - 1), (3.3) 

for every t G [0,T]. Combining (13. 3 p and (12.21) . we discover 

a 2 \\Vu Q (t)\\ 2 L2 < a 2 \\ Vu 0 || 2 2 + 2Ca||u 0 ||L2(e G:r — 1)^ 2 — C 2 a 2 (e CT — 1). 

Thus, limsup Cf _ ) . 0 +sup tg [ 0 T ] a 2 ||Vu“(t)|| 2 2 = 0, which contradicts assump¬ 
tion (13.ip . and therefore the solution u, of the 3D Euler equations, is singular 
within the interval [0, T]. □ 

3.1. Comparison with original criterion 

We show that the new blow-up criterion (13.ip is stronger than (II. ip . Since 

sup a 2 ||Vu(t )|| 2 2 > a 2 \\Vu(t)\\ 2 L2 , (3.4) 

te[o,T\ 

for any t G [0, T], we may take the limsup a _ > . 0 + of both sides to obtain 

limsup sup a 2 ||Vu(f)|| 2 2 > limsupa 2 ||Vu(f)|| 2 2 . (3.5) 

a—»0+ te[0,T] a^0+ 

The left-hand side is constant, and the right-hand side depends on t. Thus, 
limsup sup a 2 ||Vu(t )|| 2 2 > sup limsupa 2 ||V'u(t)|| 2 2 . (3.6) 

a->0+ te[0,T] te[0 ,T] a—>-0+ 

Therefore, if the right-hand side is positive, the left-hand side is positive. 
Hence, (11.2p implies (13.11) . 
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